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Abstract 

We study the equilibrium thermodynamics of a simple, confining, DSE-model of 
2-flavour QCD at finite temperature and chemical potential. The model has two 
phases: one characterised by confinement and dynamical chiral symmetry breaking; 
and the other by their absence. The phase boundary is defined by the zero of 
the vacuum-pressure difference between the confined and deconfined phases. Chiral 
symmetry restoration and deconfinement are coincident with the transition being of 
hrst order, except for p = 0, where it is second order. Nonperturbative modifications 
of the dressed-quark propagator persist into the deconfined domain and lead to a 
dispersion law modified by a dynamically-generated, momentum-dependent mass- 
scale. This entails that the Stefan-Boltzmann limit for the bulk thermodynamic 
quantities is attained only for large values of temperature and chemical potential. 
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The existence of a “quark-gluon plasma” [QGP] phase of QGD is well estab¬ 
lished. It is characterised by deconhnement, which means that quarks and 
gluons have mean free paths that are long in comparison with standard nu¬ 
clear radii, and the realisation of chiral symmetry in the Wigner mode; i.e., all 
chiral symmetry breaking effects require an explicit source: the current-quark 
mass. However, the details of the transition to this phase, and the thermody¬ 
namic parameters that characterise it: energy, entropy, pressure, etc.; remain 
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uncertain. These are the quantities that will determine whether, and under 
what conditions, the QGP can be engineered in contemporary experiments. 

The thermodynamic properties of a theory are completely determined by 
its partition function [generating functional], and this is the quantity esti¬ 
mated directly in numerical simulations of lattice-QCD actions at £nite-T. 
The Dyson-Schwinger equations [1] [DSEs] provide another nonperturbative 
means of calculating the partition function: via the evaluation of the Schwinger 
functions, which are the moments of the measure that dehnes the partition 
function. Tractable, quantitative studies of the coupled system of DSEs, exam¬ 
ples of which are the QCD gap equation and meson Bethe-Salpeter equation, 
require a truncation of the tower of integral equations; for example, via a 
choice for the quark-antiquark scattering kernel. Although it is not possible to 
judge the hdelity of a particular truncation a priori, a systematic approach [2] 
allows one to address this issue in a direct, practical and constructive manner. 
The DSEs have been applied extensively at zero temperature and chemical 
potential to the phenomenology of QCD [3], including the study of strong 
interaction contributions to weak interaction observables [4]. 

A robust, qualitative prediction of DSE studies in QCD is that, in the infrared; 
i.e., for small spacelike-g^, the propagation characteristics of the elementary 
excitations are much altered from that inferred from perturbation theory. 
For example, the dressed-gluon propagator (2-point function) is strongly en¬ 
hanced [5], which leads via the QCD gap equation to an infrared enhancement 
of the light-quark mass-function and, as a result, an infrared suppression of 
the dressed-light-quark propagator [6]. These modihcations are intimately re¬ 
lated to conhnement and dynamical chiral symmetry breaking and mean that 
new challenges arise in applying the DSEs to QCD at £nite-T and p [QCD^j. 
Therefore, as a precursor to sophisticated DSE analyses of QCD^, it is useful 
to explore simple DSE-models in which these features are manifest. 

A particularly simple and useful illustrative DSE-model of QCD [7] is specihed 
by a model dressed-gluon propagator that can be generalised to £nite-T as [8] 

Dfc) = ^ 2 ) ^ , (1) 


where (p^) = (p, D^), = 2kTiT is the boson Matsubara frequency^ and 

p is a mass-scale parameter. The infrared enhancement of the dressed-gluon 
propagator suggested by Refs. [5] is manifest in this model. As an infrared- 
dominant model that does not represent well the behaviour of D^y{p, Qk) away 
from IpI^ -|- 0, there are some model-dependent artefacts in our study. 

However, there is signihcant merit in its simplicity and, since the artefacts 

^ In our Euclidean formulation { 7 ^, 7 ;^} = 26 with 7 I = 7 ^. 
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are easily identified, the model remains useful as a means of elucidating easily 
many of the qualitative features of more sophisticated Ansatze. 

Using (1) and the “rainbow-approximation”: r“(/c,p) = for the dressed 

quark-gluon vertex; the QCD^ gap equation, or DSE for the dressed-quark 
propagator, is [2] 

S~\p,uJk) = SQ^{p,uJk) + i?7V5'(p,cUfc)7^, (2) 


where Sq^{p, Uk) = ■ p + + *h) + "rn, Uk = {2k + l)7rT is the fermion 

Matsubara frequency, p is the chemical potential and m is the current-quark 
mass. A simplicity inherent in (1) is now apparent: it allows the reduction of 
an integral equation to an algebraic equation, in whose solution many of the 
qualitative features of more sophisticated models are manifest, as will become 
clear. The solution of (2) has the general form 


S{Pk) 


1 

i7 ■ pA{pk) + ij4{(^k + ip)C{pk) + B{pk) ’ 


(3) 


where pk = (p, cufc -f ip), and (2) entails that the scalar functions introduced 
here satisfy 


= B^ + mB^ + (Apl — B^ 

2B{pk) 


m 


(2 


p^ + m^ + 4pl^ B , 


A{pk) = C{pk) = 


m + B{pk) 


(4) 

(5) 


Herein we are particularly interested in the chiral limit, m = 0. In this case (4) 
reduces to a quadratic equation for B{pk), which has two qualitatively distinct 
solutions. The “Nambu-Goldstone” solution, for which 


j Pf - 4p|, S(p|) < 4 

B{Pk) — s 

I 0, otherwise 

I 4 (^1 -h , otherwise , 


( 6 ) 

(7) 


describes a phase or mode of this model in which: 1) chiral symmetry is dy¬ 
namically broken, because one has a nonzero quark mass-function, B{pk), in 
the absence of a current-quark mass; and 2) the dressed-quarks are conhned, 
because the propagator described by these functions does not have a Lehmann 
representation. The alternative “Wigner” solution, for which 
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B{pk)=0 

(8) 

C{Pk) = 1 (^1 + ^1 +j , 

(9) 


describes a phase of the model in which chiral symmetry is not broken and 
the dressed-quarks are not conhned. 

With these two “phases”, characterised by qualitatively different, momentum- 
dependent modihcations of the quark propagator, this DSE-model of QCD^ 
can be used to explore chiral symmetry restoration and deconhnement, and 
elucidate aspects of the methodology of such studies. 


The pressure is obtained directly from the partition function, which is the 
sum of all vacuum-to-vacuum transition amplitudes. In “stationary phase” 
approximation, the partition function is given by the tree-level auxiliary-held 
effective action [9] and the pressure is: 


P|S| = ^hiZ^^ {Tl'Ln - iTl' |SS]} , 


( 10 ) 


where /? = 1/T and the self energy S(pk) = S~^(pk) — SQ^(pk). The pressure 
is a functional of S{pk)- In the absence of interactions S = 0 and (10) yields 
the free fermion partition function. 

We have neglected the gluon contribution to the pressure because, using (1), 
it is a temperature independent constant. More sophisticated Ansatze do not 
suffer this defect [8,10], which is associated with the fact that (1) does not 
represent well the ultraviolet behaviour of D^y{k) in QCD^. 

The contribution to the partition function of hadrons and hadron-like corre¬ 
lations is also neglected in (10). At the level of approximation consistent with 
(10) these terms are an additive contribution that can be estimated using the 
hadronisation techniques of Ref. [11]. After a proper normalisation of the par¬ 
tition function; i.e., subtraction of the vacuum contribution, they are the only 
contributions to the partition function in the domain of conhnement. They 
are easy to calculate and we consider them no further herein. 

In studying the phase transition one must consider the relative stability of the 
conhned and deconhned phases, which is measured by the (T, /i)-dependent 
vacuum pressure difference (or “bag constant” [12]) 


i3(T,/i)=P[ANG]-P[^w], (11) 

where S'ng means (3) obtained from (6) and Aw, (3) obtained from (8). 
B(T,p,) > 0 indicates the stability of the conhned (Nambu-Goldstone) phase 
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0.30 



Fig. 1. The phase boundary in the (T,/i)-plane obtained from (12) and (13); 
rj = 1.06 GeV [7]. The “structure” in this curve, apparent for small-T, is an arte¬ 
fact of the fact that (1) improperly represents the quark-quark interaction in the 
ultraviolet. 

and hence the phase boundary is specified by that curve in the (T, /i)-plane 
for which 


( 12 ) 

In the chiral limit, the deconfinement and chiral symmetry restoration tran¬ 
sitions in this model are coincident. The critical line is illustrated in Fig. 1, 
where explicitly 

B{T,fx)= (13) 



with; T = T/ri, n = n/rj] and i T/i^, pi = (f, Ui + ip). 


For p = 0 the transition is second order and the critical temperature is 
T° = 0.159 ? 7 , which, using the value of ?7 = 1.06 GeV obtained by htting 
the TT and p masses [7], corresponds to T° = 0.170 GeV. This is only 12% 
larger than the value obtained in Ref. [8] and the order of the transition is 
the same. However, in the present case the critical exponent is /? = 0.5, which 
differs from the result (3 ~ 0.33 obtained in Ref. [8]. This is an artefact of 
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Fig. 2. B{T,^) from (13); B{T,^) > 0 marks the confinement domain. 


the fact that ( 1 ) neglects short-range contributions to the quark-quark inter¬ 
action and, in particular, it is due to the sharp “cut-off” provided by B{pk) 
in ( 6 ). Nevertheless, the comparison with Ref. [ 8 ] illustrates that this simple 
model can provide a reasonable guide to the thermodynamic properties of 
more sophisticated DSE-models of QCD^. 

A feature of the simplicity of the model is that it allows a straightforward 
analysis of /i 7 ^ 0, where numerical simulations of lattice-QCD actions are not 
tractable [13], thereby providing valuable insight into deconhnement in this 
region. For any fi ^ 0 the transition is hrst-order, as revealed by close scrutiny 
of Fig. 2. For T = 0 the critical chemical potential is /i[] = 0.3 GeV. We also 
see from Fig. 1 that Bc{T) is insensitive to T until T 5 ^°. 

The (T,/r)-dependent vacuum pressure difference, is illustrated in 

Fig. 2. The scale is set by -B(0,0) = (0.102 = (0.109 GeV)^, which can be 

compared with the value ~ (0.145 GeV)"^ commonly used in bag-like models 
of hadrons [11]. This pressure is associated with the rearrangement of the 
conhned-quark vacuum; i.e., with the rearrangement of the “ground state”. It 
does not contribute actively to the thermodynamic pressure, describing only 
the change in the “stationary point” that dehnes the auxiliary-held effective- 
action. 

Once in the deconhnement domain, illustrated clearly in Fig. 1, the quarks 
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contribute an amount 


P[Sw]=V^‘^N,Nf 


T 


71^ 




-1 + 3? 



(14) 


to the pressure, which we renormalise to zero on the phase boundary. Just 
as in the case of free fermions, this expression is formally divergent and one 
must isolate and dehne the temperature-dependent, active contribution. This 
is made difficult by the fact that, in general, C{pi) is only known numerically, 
and hence it is not possible to evaluate P[5'w] analytically. We have developed 
a method for the numerical evaluation of (14). 

Consider the derivative of the integrand in (14): 


I i 

oo r 1 

E -i 

z=o f 


(y - P? 


+ 


{y + Pf 


(15) 


iy - pf + uf iy + /i)2 + o;, 


+ 3? 


(2C{Pi) - 1 dC{Pi] 


\ CiPiV df 


In the absence of interactions C{pi) = 1, the second term is zero and 


O OO 

/=o 


{y - pf 


{y + P? 


+ 


{y-pf+uf {y + pY + ut 


= i|f^+X(e(.))}.(16) 


where in this case e{y) = y and 


J(C) = In 


1 + exp 


C-p ' 

T 


+ In 


1 + exp 


C + 
f 


(17) 


Appropriately inserting (16) for the parenthesised term in (14), and neglecting 
T-independent terms one obtains. 


T ^ 

P[So]=v^N,Nf^ [ dyy^I{y) (18) 

71 J 
0 

=y (y+. ( 19 ) 

which is the massless, free particle pressure. 

To proceed we assume that the nontrivial momentum dependence of C{pi), 
which is manifest in all DSE-models of QCD, acts primarily to modify the 
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Momentum y 

Fig. 3. K{y, jl), which describes the nonperturbative modification of the free particle 
dispersion law, for /U = 0,0.3,0.6. By assumption, it is independent of T. 

usual massless, free particle dispersion law. We evaluate numerically the sum 
on the right-hand-side of (15) and use the form on the right-hand-side of (16) 
to £t a modified, T-independent dispersion law, e{y,jj,) = y + K{y,jl), to the 
numerical results. The existence of a K,{y, Jl) that provides a good y^-ht on 
the deconfinement domain is understood as an a posferion justihcation of this 
assumption. In our calculations, on the entire T-domain, the relative error 
between the £t and the numerical results is < 10%. 

We illustrate the calculated form of n{y,Jl) in Fig. 3; it only depends weakly 
on Jl. The form indicates that the result of the persistence of nonperturbative 
effects into the domain of deconfinement; i.e., the nontrivial momentum depen¬ 
dence of C{pi) and its slow evolution to the asymptotic value C{pi) = 1, is to 
generate a mass-scale in the massless dispersion law: fi;(0, 0) ~ 0.6 ~ 2/2°. This 
mass-scale is unrelated to the chiral-symmetry order parameter, B{0,uo + iy), 
and is a qualitatively new feature of this study. For p > 5/2° the explicit 
mass-scale introduced by the chemical potential overwhelms this dynamically 
generated scale. 

With this result (14) becomes 


P|Sw| = y AfcAf/T y dyy^I(e(y,n)), 


( 20 ) 


and the quark pressure in this DSE-model of QCD^ is 

p,(r, a) = «(!>) {P[Sw] - P[Sw]U} 


( 21 ) 








Fig. 4. The quark pressure, Pq{T,fi), normalised to the free, massless (or Ul¬ 
tra-Relativistic) result, (19). 

where V is the domain marked “Deconhned” in Fig. 1, 0{V) is a step func¬ 
tion, equal to one for (T, /i) G T>, and i^[*S'w]|gx) indicates the evaluation of this 
expression on the boundary of T>, as dehned by the intersection of a straight- 
line from the origin in the (T, /r)-plane to the argument-value. It is plotted in 
Fig. 4, which illustrates clearly that in this model the free particle (Stefan- 
Boltzmann) limit is reached at large values of T and ft. The approach to this 
limit is slow, however. For example, at T ~ 0.3 ~ 2f», or /i ~ 1.0 ~ 

( 21 ) is only f of the free particle pressure, (19). A qualitatively similar result 
is observed in numerical simulations of lattice-QCD actions at £nite-T [14]. 
This feature results from the slow approach to zero with y of illus¬ 

trated in Fig. 3, and emphasises the persistence of the momentum dependent 
modihcations of the quark propagator. 

With the dehnition and calculation of the pressure, Pq{T,y), or equivalently 
the partition function, all of the remaining bulk thermodynamic quantities 
that characterise our DSE-model of QCD^ can be calculated. As an example, 
in Fig. 5 we plot the “interaction measure”: A = e — 3P, where e is the energy 
density. This is zero for an ideal gas, hence the nomenclature: A measures the 
interaction-induced deviation from ideal gas behaviour. This figure provides 
a very clear indication of the persistence of nonperturbative effects into the 
deconhnement domain. 

We have explored the equilibrium thermodynamics of a simple, confining, 
DSE-model of QCD^, dehned by the model dressed-gluon propagator in (1), 
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Fig. 5. The “interaction measure”, A{T,n), normalised to the free, massless result 
for the pressure, (19). 

which shares many of the features found in DSE-models with more realistic 
Ansatze for the dressed-gluon propagator [8]. The simplicity of the model al¬ 
lows a straightforward analysis of the bulk thermodynamics for finite chemical 
potential as well as finite temperature. 

Deconfinement and chiral symmetry restoration are coincident and the phase 
boundary, defined by that curve in the (T, /i)-plane for which the temperature- 
and chemical-potential-dependent “bag constant”, (13), vanishes, is illustrated 
in Fig. 1. The transition is first order except for p = 0. We expect these results 
to be qualitatively reliable; a more sophisticated model leads only to a slight 
reduction in the /i = 0 critical temperature [8] and an increase in the T = 0 
critical chemical potential [10]. 

The quark pressure is given by (21) and is illustrated in Fig. 4. It is zero 
in the confining domain and approaches the Stefan-Boltzmann limit for a 
massless, free particle only for relatively large values of T and fi. The cal¬ 
culation of the pressure is complicated by the momentum dependence of the 
nonperturbatively-dressed quark propagator. This generates a nonperturba- 
tive mass-scale in the dispersion relation for the deconfined quark. Fig. 3. It 
is an essential feature of this study and is responsible for the slow attainment 
of the Stefan-Boltzmann limit. 

The persistence of nonperturbative effects into the deconfining domain is elu¬ 
cidated in the “interaction measure”, which is illustrated as a function of T 
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and /i in Fig. 5. The slow attainment of the Stefan-Boltzmann limit at finite-T 
observed in lattice simulations of QCD^ is mirrored in a slow approach to this 
limit with increasing /r. 
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